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Abstract 

We introduce an analytical solution to the one of the most familiar problems from the elementary 
quantum mechanics textbooks. The following discussion provides simple illustrations to a number 
of general concepts of quantum chaology, along with some recent developments in the field and a 
historical perspective on the subject. 
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I. INTRODUCTION 



Many aspects of the behavior of quantum systems can be understood and interpreted in 
terms of the dynamical characteristics of their classical counterparts. It is often possible 
to obtain quantitatively such important attributes of a quantum system as its spectra and 
even its wave function in terms of the appropriate classical quantities 0,0,Q|. Although the 
main objects of the classical dynamics, the dynamical trajectories, are not something that 
can be rigorously used in the context of quantum mechanics, they can often facilitate our 
understanding of quantum realm, by providing means for semiclassical interpretation. 

As an illustration of this point, let us first look at the infinite square well problem. Let 

us consider a particle confined in the infinite square well potential, 

r 0, if < X < 6, 
V{x) = \ (1) 
L oo if X < or X > L. 

Classically, the dynamics of such particle is as simple as it can possibly be - the particle 

simply bounces periodically between the two walls. Geometrically its trajectory is a closed 

loop, which the particle traverses over and over again. The total length of the loop is twice 

the width of the well, C = 2L, and the period of motion is T = 2L/v, where v is the speed 

of the particle. As a result of the wall reflections, the momentum of the particle changes its 

sign, but its magnitude never changes. 

Quantum mechanically, this system is nearly as simple to describe. The wave function 

■^(x) in this case is a combination of two plane waves, 

ij{x) = Ae'^'' + Se-*^^ (2) 

where k = \fE/h is the particle's momentum, which must satisfy the boundary conditions 
'?/'(0) = ip{L) = 0. This leads to the quantization condition 

sin(A;L) = 0, (3) 

and hence the quantum spectrum /c„ of this problem is given by 

Lkn = nn. (4) 

Note, that the left hand side of this equation coincides with the classical action integral, 

S = kdx = Lk, (5) 
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taken along the trajectory 7 that connects the two turning points, and hence the condition 
(13)) imphes that the action takes only discrete values, 5*^ = irn. Expressing the de Broglie 
wavelength A in terms of the momentum, A = 2TT/k, the relationship (0]) can be cast into 
the form £/A„ = n, which has a simple geometrical meaning. Apparently the quantization 
condition (@]) implies that the wave (j21) must fit geometrically n times onto the (only) classical 
periodic orbit that has the length C As the energy of the particle increases, its wavelength 
will become smaller and smaller, however the geometrical condition £/A„ = n always holds. 

Such simple combination of physical and geometrical ideas were used by Bohr and his 
school in around 1914 to provide the first explanations to one of the most striking features of 
the quantum mechanical systems - the discrete nature of their energy spectra. According to 
their views, the discreteness of the quantum spectra turns out to be essentially a consequence 
of the geometrical consistency of the wave mechanics. Specifically, such "wave-geometrical" 
approach proved very successful in early attempts to explain the experimentally observed 
emission-absorption spectrum of the Hydrogen atom. In fact, the results obtained in this 
way were exact, and that gave serious reasons to believe that these ideas were adequate to 
describe the quantum physics of the subatomic world. 



However, the later attempts of Bohr, Sommerfeld, van Vleck, Born and others 
to continue the success of these ideas on more complicated atoms, failed completely. 
It was apparently impossible to explain even approximately the experimentally observed 
spectrum of the Helium atom - the next simplest atom after Hydrogen in the periodic table 
of elements and hence the next best candidate for a successful treatment by means of such 
wave-geometrical quantization. Van Vleck wrote in 1922 j^: 

''The conventional quantum theory of atomic structure does not appear able to account 
for the properties of even such a simple element as helium, and to escape from this dilemma 
some radical modification in the ordinary conceptions of quantum theory or of the electron 
may he necessary" . 

What was be the nature of the difficulties that required such "radical" and "conceptual" 
modifications? In 1925 Max Born wrote ^: 

"...the systematic application of the principles of the quantum theory... gives results 
in agreement with experiment only in those cases where the motion of a single electron is 
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considered; it fails even in the treatment of the motion of the two electrons in the helium 
atom. 

This is not surprising, for the principles used are not really consistent... A complete 
systematic transformation of the classical mechanics into a discontinuous mechanics is the 
goal towards which the quantum theory strives. " 

Seemingly one makes a very simple and natural move by trying to go from the exactly 
solved Hydrogen atom to the Helium, by adding just one more particle to the two body 
nucleus-electron system of the Hydrogen. However, from the point of view of the contempo- 
rary classical mechanics, this modest generalization turns an integrable two-body system into 
a non-integrable three body system. Hence, in order to impose the "geometrical consistency" 
on the quantum waves that could describe the Helium atom in Bohr's approach, one would 
have to deal with the overwhelming geometrical complexity of its classical phase space. It 
means that all the familiar "niceties" of the generic chaotic systems, such as the exponential 
proliferation of the periodic orbits associated with extreme complexity of their shapes, would 
have to be taken into account - something that could hardly had been done at Bohr's time. 
Indeed, at the beginning of the last century chaos theory was just being developed in the 
works of Poincare, Lyapunov, Hadamard, Birkhoff and a few others Q^laSll^l- Although 
the importance of the classical dynamical behavior for successful quantization was realized 



by certain researchers ll| , it took about 60 years before the first semiclassical quantization 



procedure for nonintegrable systems was outlined. As for the Helium atom, a semiclassical 



quantization scheme for it was proposed in 1992 
quantum chaos theory P|. 
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using some recent developments of the 



II. BACK TO ID POTENTIAL WELLS 



Surprisingly, the difficulties of the semiclassical quantization of the classically noninte- 
grable Helium atom system can be illustrated by means of elementary quantum mechanics. 
Below we shall consider a simple modification of a completely transparent system (Q) that 
leads to a transition from classical integrability to a non-integrability and clarifies the cor- 
responding outburst of complexity on quantum level. 

Let us add a small step at the bottom of the potential (jT]) and consider a point particle 
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FIG. 1: A simple modification of tlie infinite square well model - a square well with a potential 
step, along with the two apparent classical trajectories 71 and 71. 

moving in the potential 

{0, if < X < 6, 
V, iib<x<L, (6) 
00 if X < or X > L, 

shown in Figure 1. Seemingly the classical picture doesn't change much. It appears that 
if the energy of the particle is higher than V, the particle will oscillate between the points 
X = and x = L just as before, and if its energy is below the potential step hight, it 
will oscillate between x = and x = b. However, the reality turns out to be much more 
complicated than that. 

Let us first give the quantum mechanical description of the problem. Suppose that the 
energy E of the particle is above the potential height V. The momentum of the particle in 
the region < x < b is k = ^/E, and in the region 6<x<litisK = \/E — V, and its tjj 
function consists of two parts, 

{sm(kx), if < X < 6, 

- - 
A sin [k{L — x)] , if b < X < L, 

which match continuously at x = 6. From this continuity requirement one can determine 
the spectral equation for the spectrum of the problem, 

sin(Lifc + L2K) = r sin [{Lik - L2K)] , E > V, (8) 



where Li = b and L2 = L — b are correspondingly the lengths of the left and the right sides 
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FIG. 2: A quantum particle in a square well with a step. 



of the well (jHI) and r is the reflection coefficient 

k — K 



(9) 



k + K 

The E < V case can be treated similarly and basically amounts to substituting n ik. 



X = \/V — E, in the expressions (jH} and After extracting the complex phase, the 
spectral equation becomes 

sin Life = y ^ K J E <V, (10) 

where n = 1, 2, is the root index. 

There is an important analogy between (jHJ, ()10j) and Q - the arguments of the sine func- 
tions in the left hand sides of these equations are the action lengths S{E) of the classically 
available regions in the well ^ for E > V and E < V correspondingly, 

r ( kLi + kL2, if E >V, 

S(k) = kdx={ (11) 

Uli he <V, 

which is a continuous and monotonically increasing function of the energy. 

On the other hand, it is interesting that unlike Q, the equations (jH)) for E > V and (fTUj) 
for E < V are transcendental equations - that is substantially more difficult to solve. In 
fact, in absence of any analytical ways of solving the equation (jH]) and (|Tn|l explicitly, all the 
standard textbooks, e.g. [l^, JJ, use graphical or numerical methods to approximate 

its roots, kn. Notably, there are also no simple "semi classical" interpretations similar to Q 
for the results of such approximations, which, after all, does not appear important for the 




FIG. 3: Typical appearance of the graphical methods currently used by the Quantum Mechanics 
textbooks. The roots are found as the intersection points of the graphs representing the right and 
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eft hand sides of a suitable spectral equation, marked by the vertical dashed lines. For details 
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purpose of obtaining the numerical values for the roots of the spectral equations (jHI) and 

This is very similar to what happened to the Helium atom quantization problem - the "old 
quantum theory" treatment became unnecessary once the formalism of the "new" quantum 
mechanics of Schrodinger and Heisenberg was established and could be used. The problem 
of obtaining the energy spectrum of the Helium (or in principle any other atom) was reduced 
to a technical problem of diagonalizing the Hamiltonian matrix Hnm by means of all sorts 
of numerical techniques and approximations. 

Even without making any historical references, one can notice the obvious contrast in the 
level of complexity between the spectral equations (jH)), (fTUI) and Q. What is the physical 
reason for it? Below we shall argue that the complexity of the spectral equations (jH)), ()10|) 
has a deep physical meaning and can be understood from analyzing the classical motion of 
a particle in the potential well 
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III. THE CLASSICAL LIMIT 



To start the discussion of the classical dynamics of a point particle in the potential ©, 
let us notice that the reflection coefficient Q in (jH)) (i? > does not depend on the Plank's 
constant h, 

1 - J\ - VIE , , 

r = ^ , ' . (12) 

1 + V^l - VIE 

Due to this circumstance this quantity (although obtained by purely quantum mechanical 
means is in fact quite classical, because one does not need to refer to any quantum 
mechanical concepts to evaluate r. 



This curious fact can be easily interpreted and understood from general principles 17|, 3| ■ 



Just as any other wave, the probability amplitude (j21) gets reflected or diffracted when it 
encounters inhomogeneities on its way. It is important that the scale of these inhomogeneities 
should not be too big compared to the wavelength of the wave, otherwise the wave will 
"adjust" to the smooth changes of the properties of the media. It is fair to say that in 
order to induce reflections, the obstacle should appear somewhat "abruptly" in front of the 
wave, at a scale d smaller than the characteristic wavelength, d \. In our case, since 
the potential step © is defined to be absolutely sharp, changing discontinuously at x = 6 
from to V ^ the quantum-mechanical wave will be always reflecting from its boundary at 
X = b with the reflection probability r^, no matter how small its wavelength is, even at the 
classical limit A = 0. 

In other words, although intuitively one would expect the particle moving with the energy 
E > V only to change its speed after passing over the point x = b, there actually exists a 
possibility of classical (the so called non- Newtonian, [itI ll^) backward reflections from the 
sharp potential barrier edge. 

This curious fact emphasizes an interesting aspect of the connection between the classical 
and quantum mechanics, known as the "Correspondence Principle" , which was first invoked 
by Niels Bohr in around 1923. This fundamental principle states that classical mechanics 
can be understood as a limiting case of quantum mechanics in the so called "classical limit" , 
i.e. in case when the motions are characterized by the actions much larger than the value 
of the Plank's constant h. 

At Bohr's time, the fulfillment of such quantum-classical correspondence was viewed as 
a natural way to validate meaningfulness and physical consistency of the quantized analogs 
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of familiar classical systems, such as atoms. Hence the correspondence principle is often un- 
derstood as a naive requirement for the quantum system to reproduce the expected classical 
behavior in the limit h ^ 0, whereas in fact, there is no such requirement. The appearance 
of Non-Newtonian scattering events is a curious example of a situation when quantum me- 
chanics elucidates a certain implicit aspect of the corresponding classical dynamics, bringing 
up details that could have been easily overlooked. 

Similar considerations can prove that non-Newtonian scattering phenomena can happen 
for every potential with sharp edges, as some kind of a reminder of quantum-mechanical 
legacy of classical mechanics. The fact that we never observe such events in every day life 
implies that in reality there are no potentials changes, sharp on the scale of the quantum 
wavelengths of the macroscopic objects. 

IV. NON-NEWTONIAN CHAOS 

Because of the possibility of such classical non-Newtonian reflections, the classical dynam- 
ics of a particle in the potential © is far from trivial. Every time the particle approaches 
the boundary x = b between the two regions in ©, it can be reflected from it with the 
probability and transmitted through it with the probability t"^ = 1 — r^. As a result, 
instead of a couple of back and forth oscillations which one would naively expect from a 
particle in the potential (jHl), the actual trajectories of such particle are far more complex. 
The particle can start, say, in the left side L of the well, move to the boundary x = b, reflect 
back from it, reflect back from the rigid wall at x = 0, do this several times, then transmit 
eventually to the right side R of the well, where it can also perform a number of oscillations 
before returning to the left side of the well, etc. Any thinkable sequence of oscillations in 
the right and in the left sides of the potential well (jH)) represents a possible trajectory of the 
particle. It should also be emphasized, that at every reflection or transmission (scattering) 
event, the particle completely looses its memory about the previous stage of its motion. 
Given the current position and the momentum of the particle, neither its previous evolution 
nor its state of motion after the next collision can be reconstructed. Hence, instead of the 
deterministic evolution we generated a fairly complicated stochastic dynamical process by 
considering a seemingly simple potential (jUj). 

Obviously, there can be classical orbits of arbitrary length, and the bigger is the allowed 
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FIG. 4: A storm in a coffee cup - non Newtonian chaos in a one dimensional step potential. The 
lines represent examples of non-Newtonian orbits in the potential For illustration purposes 
different trajectories are shown at different energies, as well as different sections of same trajectories 
are shifted slightly up and down the energy axis, so that their parts are visible. 

length (or the period) of the trajectory, the larger are the numbers of nontrivial orbits that 
come into picture. Is there a way to enumerate all this variety the orbits? It turns out 
that the description of the general behavior of the orbits in this system can be conveniently 
formalized. Indeed, every orbit can be described by a two-letter code, which would simply 
tell us in what sequence the orbit swings through the left (L) or the right (i?) sides of the 
well, as shown in (Fig. 4). Periodic orbits are obviously represented by periodic sequences of 
symbols. Conversely, any periodic sequence of the symbols (for example, RLRLRLRLRL... 
or just {RL) for short), unambiguously describes a certain periodic orbit (note however, 
that two sequences which can be obtained from one another by a cyclical permutation of 
symbols, correspond to the same periodic orbit). If the orbit swings rii times in the left side 
of the well and n2 in the right side, its action length is equal 

Sorh = "^riiSi + 2n2S2, (13) 

where Si = Lik and 5*2 = L2K. 

The number of the periodic orbits increases indefinitely. Their shapes (and correspond- 
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FIG. 5: A couple of examples of a non-Newtonian orbits in the potential © marked by their 
R-L binary codes. Every left swing of the orbit contributes an "L" to its code and every right one 
contributes an "R". 

ingly their L—R binary codes) become more and more complicated, so the classical mechanics 
in the potential turns out to be surprisingly rich. In fact, the number of the geometrically 
different orbit shapes (or the number of the prime periodic orbits, the ones that never retrace 
themselves) that include up to m scattering events grows exponentially as 

p0.7m 

N ^ . (14) 

m 

Such behavior closely resembles the periodic orbit proliferation scenario that takes place 
in classically nonintegrable chaotic systems. Hence a particle in the potential © can be 
viewed as a simple model of a classically chaotic system. This L — R code representation 
of the orbits is a simple example of symbolic dynamics over a partition of the phase space 

QQ. 

The difference between this and the simple back and forth motion in the square well 
potential is overwhelming, and provides a simple illustration to the difference between the 
dynamical behavior of the integrable and the nonintegrable systems, such as e.g. Hydrogen 
and Helium. The outburst of complexity that results after adding an extra potential step 
at the bottom of the potential is similar to the situation with the Helium atom. An 
extra electron that is added to the two-body system of the nucleus and the electron in the 
Hydrogen atom totally destroys its integrability. 

Of course the more complex is the orbit shape, the smaller is the probability that this 
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FIG. 6: Exponential proliferation of the orbits in the potential Q. The orbits that undergo 1, 2, 
3 and 4 scatterings marked by their R-L binary codes. 

particular orbit will be realized. This "complexity selection" can be easily understood from 
the point of view of the quantum mechanics, where the propagation of the particle is de- 
scribed by the normalized wave It is clear that every reflection or transmission reduces 
its initial amplitude by the amount of the reflection or transmission coefficient. Hence, it 
is clear that if a certain prime periodic orbit p (i.e. the one that can not be considered a 
repetition of a shorter orbit) reflects a{p) times from the either side of the potential barrier 
at X = b, and transmits r(p) times through it, its initial amplitude will decrease 

times Here the factor (— keeps track of the sign changes due to the wall reflections 
and the right reflections from the boundary. If the orbit p is traced over itself u times, then 
the corresponding amplitude will be A^. So the quantity (fT^ is a certain quantum (and, 
since r does not depend on h, also classical) "weight" of the orbit. 
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V. CLASSICAL DYNAMICS AND QUANTUM SPECTRUM 



In the case of the infinite square well potential the orbit bouncing back and forth between 
the walls of (P) completely exhausts all the classical dynamical possibilities of the particle 
in the infinite square well. This orbit also defines the quantum energy spectrum (0]) of 
the particle via the relationship Sn = 7m. This is analogous to the general quantization 
procedure of the so called integrable systems, i.e. the ones that have as many dynamically 
conserved quantities (e.g. angular momentum, etc.) as degrees of freedom jll|, |25|. These 
conserved quantities (the action integrals) /«, i = l,...,d, are quantized according to 

li = h{ni + fii), (16) 

where n^s are natural numbers, = 1,2, and /ij is a certain geometrical constant jioj ]. 
The energy of these systems is a certain function of the Jj's, E = E{Ii, ...In), and so the 
quantization rules ()16|) provide immediately the (semiclassical) energy spectrum. For exam- 
ple, for a particle moving in the infinite potential well (^, the magnitude of the momentum 
and hence the action, / = 2|p|L, is conserved. It is then quantized as \p\ = hnn/L, and 
yields the quantum energy eigenvalues E = (/i7m)^/2mL^. 

According to this scheme, the "quantum numbers" rij are naturally related to the classical 
integrals of motion Jj. 

On the other hand, the classical behavior of the particle in the potential © is of manifestly 
nonintegrable type. The variety of possible dynamical orbits in the potential (jH)) is much 
richer than in (^, and we were able to describe all of them. What information is hidden 
in this variety of the orbits? Is it possible, after all, to use this knowledge to obtain the 
quantum energy spectrum? One should remember that most of these orbits appeared as 
some leftovers of purely quantum effects - the over-barrier reflections and under-barrier 
tunnelings. Do they carry through the complete quantum mechanical information about 
the particle? 

As pointed out above, in the integrable systems there exists a convenient handle - the 
action integrals, which can be ascribed certain discrete values via the EBK semiclassical 
quantization procedure above 11, 2^. What if the integrals do not exist? What is there in 
a nonintegrable system that can be used in order to quantize the system semiclassically? 

The flrst answer to this question came in early 70s, when Gutzwiller [2^ showed that 
a "handle" for semiclassical quantization of the systems without the sufficient number of 
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the integrals of motion is provided by the so-called density of states - a functional that 
corresponds a 6 -peak to every energy level En of the quantum spectrum, 

oo 

p{E)^Y.^{E-En). (17) 

n=l 

Gutzwiller provided a semiclassical expansion for p{E) (the so-called Gutzwiller trace for- 
mula) in terms of the classical quantities related to the periodic orbits, 

p{E) ^ p{E) + - Re 5^ T, 5^ B; e-^^(^). (18) 

p u=l 

Here Tp, Sp and Bp are correspondingly the period, the action and a certain weight factor 
(see below) of the prime periodic orbit p (all classical quantities), and u is the number of 
times the orbit p repeats itself. The first term p{E) represents the non-oscillating part of 
the density of states. 

The gist Gutzwiller 's trace formula is an (almost miraculous) interference effect, produced 
by the infinity of oscillating terms (one per each periodic orbit) in the sum (fTS|) . The 
statement made by (fTH|l is that if the energy E happens to coincide with a quantum energy 
level En, then all the terms in the sum (jl8|) will interfere constructively and produce a peak, 
whereas if there are no quantum EnS there, they will interfere destructively and yield 0. 

This is a truly remarkable connection between the classical and the quantum properties 
of a system. After all, from a formal mathematical perspective, classical characteristics 
of a system should a priori describe only its ^ = limit, whereas according to the trace 
formula one can use the classical properties of the system to extract the information about 
its quantum properties for h 0. 

Gutzwiller's formula can also be applied to the ID system of a particle in the step 



potential ©. Moreover, it can be shown 26| that for the system © the expansion ()18|) 
is exact, although usually it provides only an approximate (with semiclassical accuracy) 
representation of p{E). Another important characteristic of the spectrum is the "spectral 
staircase" N{E), 

N{E)= I p{E')dE' = y2e{E-En), (19) 
where 9(x) is Heaviside's theta function, 

f 1, ifa;>0, 

e(x)= (20) 
I li X <0, 
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which gives the number of the energy levels on the interval between and E. N{E) can 
also be expanded into a periodic orbit series, 

1 °° A'^ 

N{E) = N{E) + -lmY,Y. — (21) 

p v=l 

What makes the expansion ^Z\\ji particularly convenient for our use, is that the weight factors 
Ap in it are explicitly given by pSj). The first term of N{E) (its non-oscillating part, also 
known as Weyl's average) is given by 

N{E) = -S{E) - 70, (22) 

71 

where S{E) is the action length of the classically available region of the well, and 70 is 
a small correction term (almost a constant, see below). 

The complexity of our task of describing the spectrum in ® can already be appreciated 
from the expansions (fTHj) and ()21|) : every non-Newtonian orbit that exists in © explicitly 
contributes to these expansions. 

In addition to the huge number of the periodic orbits, there are many other subtle dynam- 
ical effects that contribute to the complexity of the expansions (jl8|) and (j21|) . For instance, it 
was mentioned above that as the energy of the particle passes from above to a level below 
V, the particle becomes classically unable to penetrate into the right section of the well, 
b < X < L. However, there exists the possibility of quantum tunneling into the region R. 
This is manifested by the corresponding momentum becoming imaginary, n ix. Hence, 
the contributions to the orbit actions (fT!^ due to the "tunneling" through the step parts are 
also imaginary. Sorb = 25*1^1 + ^2^25*2 (such orbits were called "ghost orbits" in j|2l|). 

That implies that as the energy of the particle changes, the physical characteristics of the 
orbits (and therefore of the expansion terms in ()18|) ) may change (phenomenon known as 



the "phase space metamorphosis" 



22|). However, the formal structure of the Gutzwiller's 



formula (at least in the case at hand) remains the same, and so the expansions p8|l and 
(j2Il) can be used, with due care, both for > and for E < V. 



VI. OBTAINING THE SPECTRUM 



Since the right hand side of the expansion ()18p can be obtained from considering the 
classical motion of the particle, the expansion |TH|) provides a clear connection between the 
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dynamical characteristics of the system and its spectrum in quantum regime. The energy 
levels of the quantum system are obtained according to (fTH|) as the poles (delta spikes) 
produced by the periodic orbit sum 3, [2^. It is important however, that by using this 
approach, one can not tell when these poles will appear prior to performing the summation 
of the periodic orbit series (jl8p . Without having any extra information about the system, 
the only general strategy for obtaining -E^'s is to scan the energy axis by summing the series 
fll8|) for every value of E to find out at what energies the sum produces a 5-peak. 

This illustrates the fact that p{E) is a 5^/060/ characteristic of the spectrum. Both p{E) and 
its expansion ()18j) describe the whole spectrum at once, rather than the specific individual 
energy levels, in contrast with the case of the integrable systems, where every action integral 
Jj is quantized directly and separately from the rest of the degrees of freedom. 

Actually, it is possible to extract the information about the individual energy levels 
En out of p{k) without such tedious "energy axis scanning" - if only one knows (at least 
approximately) where to look. Indeed, if for instance we would happen to know that some 
energy level En is the only level that lies between two points E^ and E' (and hence p{E) 
has only one 6-peak between these points), then we could evaluate the integral 

rE+ 

En= Ep{E)dE, (23) 

J En 

and obtain the value of En- Since we know the expansion of p{E) in terms of the periodic 
orbits (the right hand side of (fTHj) ). this integration (at least in principle) can be performed. 

If our knowledge of a certain system would be so complete that we would be able to 
separate every energy level En from its neighbors by two "separators" ^+ and E~ , En > 
En > E~, then we could use to obtain an exact numerical value of every energy level 
in the spectrum of our problem out of p{E). Apparently, for the purpose of separating one 
level from another, it is enough to assume that En-i = E', so basically we are looking for 
a sequence of points En that interweaves the sequence of the energy levels and separates 



one energy level from another 
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In other words, in order to be able to obtain 



quantitatively any energy level En out of p{E), we need to establish a partition of the E-axis 



into the intervals 



En-l, En 



each one of which contains exactly one energy level. 
Note, that even if that could be achieved, one would still have just an algorithmic recipe 
for evaluating EnS rather than a formula of the type En = ■■■ ■ In order to get such a 
formula, we would need to find a global function E{n) that depends explicitly on n and 
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produces all the separating points in their natural sequence as a function of their index, 

En = E{n). (24) 



Once such a functional dependence of the separators on their index is established, the integral 
(I2H|) would turn the index n into a quantum number and give us a complete solution to the 
spectral problem in the form En = ■■■ ■ 

The problem is that usually this is not an easy task. In order to establish a partition 



En-l, En 



one needs to have some extra 



of the energy axis into the separating intervals 
information about the behavior of the spectrum. Luckily, such information can be obtained 
for the particle in the potential (jUj). 



VII. SPECTRAL EQUATION 

In order to get this information, let us examine closely the spectral equations (jHI) and 
(ll(J|) . Since the quantity that we are after, the classical action S{E) defined by is in the 
arguments of the sines in the left hand sides of the equations (jH)) and (fTn|) . we can formally 
invert them to obtain 

/ 1\ f (-l)"(arcsin(rsin(5i-52) + f)) , ifE>V, 

S = n { n — < (25) 

V 2; \arctan(x(l + e-2^2)/^(l_e-252)) ^ ifE<V, 

where n = 1,2,.... This form of the spectral equation is particularly important for two 
reasons. First, we are getting an index n which (one would assume) numbers the solutions 
Sn - the discrete quantum values of action. 

One can notice an interesting resemblance between the way the formulae ()25|) begin, 
= TT — 1) + and the formula (fTBj) . One can speculate therefore that the first term, 
TT (n — I), expresses the regular part of the spectrum, and the second terms, which actually 
make the equation ()25j) transcendental, are introducing the irregularities into the spectrum, 
which are due to the classical non-integrability of the potential well (jHI). 

This in fact turns out to be a very important observation. Is the spectrum defined by 
fl25|) more "regular" or "irregular"? Which one of these two terms contributes more to the 
solution, the regular vr (ri — 1) part or the irregular second term? Surprisingly, it turns out 
that in terms of the magnitudes the regular part wins. Indeed, it is well known that the 
inverse trigonometric functions are bounded, — | < arcsin(a;), arctan(x) < |. On the other 
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hand, when the index n in the equation ()25p changes by 1 as we go from one level to another, 
the "regular term" increases by vr, which is generically (that is almost always) more than 
what the arcsine or arctangent can provide. 

Therefore, it is the regular first term — in (j25|) that contributes the most to the 
solution. If the step at the bottom of the well would disappear, so would the "transcendental" 
terms in ()25|). and a pure regular (periodic) spectrum Q for S would be recovered. Hence 
the second term in (pSjl can be regarded as the one responsible for the chaos induced spectral 
"fluctuations" . 

The second important consequence of writing the equations (jH)) and (jlOj) in the inverted 
form (j2SI), is that it allows one to obtain the set of separating points discussed in the previous 
section, and hence to implement the solution to the spectral problem for the potential (jH)). 

Indeed, since the irregular terms in ()25|) are different from 0, the values S'„ = 7r(r?, — 1/2) 
themselves are never the solutions to ()25|). On the other hand, the difference between Sn 
and the roots Sn-, which is due to the second term, is smaller than vr, and as a result the 
roots Sn of pKjl will be locked inside of the intervals 

vr (ri- i ) < S'„ < vr (n+ - ) . (26) 



2 J V 2^ 

Thus, the first discrete quantum action value will be locked between 7r/2 and 37r/2, the 
second one between 37r/2 and 57r/2, and so on. In other words, we run immediately into the 
separating points for the quantum values of action, 

Sn = T^{n- 1/2), n=l,2,..., (27) 

or the action separators. These separators can be used to extract the energy separators En 
from S{En) = T^i^ — 1/2), which then can be used to obtain the energy spectrum via (j^^ . 
Alternatively, rather than solving the equation S{En) = vr(n — 1/2) for the -En's in order 
to use them in the formula (PHj) . one could introduce a "density of the action states", p{S), 
defined as p{S)dS = p{E)dE, 

°° dF 

p(^) = ^<5(5-5„)^, (28) 

n=l 

to obtain the actual quantum levels of action first and then to extract the quantum energy 
levels from them. This would completely parallel the p{E) approach outlined in the previous 
section. A simple change of variables will produce the periodic orbit expansion for p{S). 
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FIG. 7: The energy separators En (their values on the energy scale are shown by solid lines) define 
the upper and the lower bounds for the quantum energy levels shown by the dashed lines. 

Hence it is possible, in accordance with the ideas outlined above, to obtain the discrete 
quantum action values Sn via 

r.7r(n+l/2) 7r(n+l/2) /•7r(n+l/2) 

Sp{S)dS = SN{S) 

'7r{n-l/2) 

where we used the identity 



Mn+l/2) 7r(n+l/2) Mn+l/2) 

Sn= Sp{S)dS = SN{S) - N{S)dS, (29) 

J-K{n--l/2) 7r{n-l/2) Jn{n-l/2) 



piS) = (30) 

Using the relationship N{Sn) = n (there are exactly n roots Sn below the nth separator) 
and the expansion (pT|l . we get 

Mn+l/2) _ . rnin+1/2) J^u 

Sn = 27171 / N{S)dS - - Im V / ^e'^^^dS, (31) 

2 A(n-l/2) TT A{n-l/2) ^ 

where the weight factors Ap are given by (fT^ . 

Note, that Q is exactly the first term of the expansion (jHT|) . Hence the oscillatory terms 
in ()31|) . with their amplitudes proportional to the non- Newtonian reflection amplitudes, are 
indeed due to the "nonintegrability" of 

Since all of the quantities on the right-hand side of ()3H1 are known, this formula provides 
an explicit representation of the discrete quantum values of the action functional in terms 
of the geometric and dynamical characteristics of the potential (jH} • Formula allows the 
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computation of the action corresponding to every quantum level individually, explicitly and 
exactly in terms of the classical parameters, indexed by its "quantum number" n. 

It should be mentioned however, that the quantum number n in (0]) and in (j31|) have 
rather different origins. Formula (j3H) does not have such a direct geometrical interpretation 
as (0]). While in (0]) the number n is the number of waves that fit on the trajectory, in the 
case of the step potential © n is the index of the cell that contains the corresponding 
action value Sn- 

So once again, we arrive at defining the quantum energy spectrum via a discrete set of 
allowed values of the action functional, Sn = S{En). This relationship (which is yet to be 
presented explicitly) can be viewed as a direct generalization of (@)). 

Of course, for using the periodic orbit expansions in practice, one needs to know how to 
truncate the series ()18|) or ()31|) to obtain finite order approximations to Sn- Since usually 
the periodic orbit series are not absolutely convergent, the order in which the expansion 
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I, that 



terms are incorporated into the sum ()3Tjl is important. It turns out 
the correct way to obtain mth approximation to the exact value of Sn is to include into the 
sum (jnH) all the trajectories that reach the point x = b m times or less. 



VIII. OVERVIEW AND AN EXAMPLE 



Let us summarize the steps for obtaining the mth correction to Sn using the semiclassical 
periodic orbit expansion technique. 

1. Write down all the m-letter sequences (words) Wm- For example, there are 2 words, 
(L) and (R), for m = 1, 4 words (LL), (LR), (RL) and (RR) for m = 2, 8 words 
(LLL), (LLR), (LRL), (RLL), (RLR), (RRL), (LRR) and (RRR) for m = 3 and so 
on. (To do the same thing geometrically - draw all possible periodic orbits similar to 
the ones shown in Fig. 5 that include exactly m transmissions and refiections). 

2. Find which w's are cyclic permutations of one another - all these sequences represent 
the same orbit, so pick one and discard its replicas. From the above examples, both 
m = 1 orbits will remain, for m = 2 we may keep (LL), (LR) and (RR), for m = 3 
we keep (LLL), (LLR), (RLR) and (RRR). (Geometrically: find out which loops 
are representing the same sequence of left and right swings. Pick one and discard its 
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replicas) . 

3. Count the number ul of L's and ur of -R's in w and find the action of the orbit w 
according to 

= 2nRSi + 2nLS2. (32) 
Remember that ii E < V then 5*2 = iy/V — EL2. 

4. Assuming that the first symbol in w cyclically follows the last one, find out which of 
the remaining sequences w are prime sequences and which ones are repetitions of a 
shorter code. Find the repetition number for each orbit (use z/^ = 1 if the orbit is 
prime). For example, for m = 2 the orbit = 1, z/^jr = ull = 2, and for m = 3 
^LLL = ^RRR = 3, I'LLR = ^RLR = 1- (Geometrically, find out how many times each 
orbit traverses over itself). 

5. Scan each word Wm (again, the first symbol of w follows the last one) and write down 
the weight according to the substitutions 

LR ^t, RL^ t, LL ^ r, RR ^ -r, (33) 

plus each wall reflection contributes a sign change. For example, for m = 1, = —r, 
Ar = r, for m = 2 Alr = t"^, Arr = All = r"^, and for m = 3 we have Alll = —r^, 
Arrr = r^, Allr = —rt"^, Ajilr = ft^- (Geometrically: assign a factor t to every 
transmission, a factor —1 to every wall reflection, a factor r to every left side reflection 
and a factor — r to every right side reflection). 

6. For every orbit, find out the corresponding term in the expansion (j!?T| . 

/ ^e'^-(^)t/^, (34) 

Jn{n~l/2) 

and integrate. 

If all the words up the length m are considered, the result will produce mth approximation 
to the exact discrete sequence of actions. 

All the non-Newtonian orbits described in the previous sections contribute to the expan- 
sion of each eigenvalue S'„. One can now appreciate the complexity of solving analytically 
the equations ^ and (fTUj) compared to solving the simple equation Q, by comparing the 
complexities of the corresponding classical dynamics. 
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FIG. 8: The orbits required for evaluating the first three periodic orbit expansion approximations 
to the Sn- Since the evaluation of the low levels which correspond to the energy E < V will require 
ghost orbits, a few ghost orbits are also illustrated. The tunneling parts of the ghost orbits are 
shown with dashed lines. 

A. An example 

As an example, let us compute the first three approximations to a a few momentum 
eigenvalues, using the orbits that include up to 3 scattering events, outlined in the summary 
and shown in Fig. 6. We have according to (|31|) that the first three corrections will be: 

^7r(n+l/2) 

Sn = 27m / N{S)dS 

/.7r(n+l/2) 

- Im / (re*2Sifc _ ^^2S,k^^ 

/.7r(n+l/2) /. . \ 

- Im / ( -r^e^^^^ + ^2^^251+^25, ^ £^4.5.^2 \ 

Jn{n~l/2) \2 2 / 

- Im r^"""'^'^ (^r'e'^'^ + tVe^(25.+450 _ ^2^^.(251+45.) _ 1^3^6.5^ ^35) 

In order to use the formula ()35|) . one needs to specify the parameters of the potential well 
(jSJ. In particular, it is a specific value of V that determines at what energy the periodic 
orbits will acquire tunneling parts. For simplicity, let us consider the case of symmetrical 
well, when Li = L2 = 1. Then the action 5* defined by (fTT| will be 

k + K, if E >V, 
k ifE<V, 



S{E) 



(36) 
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in terms of which 

r = e^^'^^o, 70 = ^ arctan , ^2 ^ (1 - e^^'^To) if E <V. 

The phase 70 (which also appears in the Weyl's average (|22p) is 70 = 1/2 for E' > V^. 

Next, let us pick a certain value for the potential step hight, say V = 25. This will 
set the Ecrit = 25 and kcru = 5 to be the "critical values" of energy and the momentum 
- that is if we will be looking for an energy level below E = V = 25, then the tunneling 
effects (ghost orbits, the restructuring of the periodic orbit sum (jl8|) ) have to be taken into 
account. Correspondingly, the critical value for the action will be Scru = kcruLi = 5, and so 
one should use different parts of formulae (jHUj) to extract the En or A;„ out of Sn, depending 
on whether a particular s„ is smaller or greater than 5. 

After the height of the potential step has been chosen, certain general statements about 
the behavior of the quantum Sn levels can be made. For the specific case V = 25, the first 
and the second separators, 5*1 = | and 6*2 = will be smaller than Scrit, and so the level 
contained between them will correspond to the energy Ei < V. The third action separator 
is bigger than kcrit, S3 = 57r/2 > 5, so we will be in a better position to judge about the 
location of the second level ^2 < 5*2 < 5*3 after computing the corrections (jHSj). All the other 
levels Sn, for n > 2, have energies higher than V. 

So let us find the discrete quantum action values 5*1, 5*2 and (for example) 5*17, and the 
corresponding quantum momentum eigenvalues, ki, k2 and kn. Using the expressions (jHSI), 
(jSni), and integrating from tt/2 to 37r/2, from 37r/2 to 57r/2 and from 337r/2 to 357r/2 
(using some numerical integration software is highly recommended) we find the values shown 
in Table 1. 

Note that in the second case the integration in ()35p should be split in two parts, 

p5tt/2 p5 p5tt/2 

/ N{S)dS= / N^^<^\S)dS+ / N^^>^\S)dS, (38) 

in order to follow the changes in the structure of the expansion terms described by (jHU)) 
and (jHTj) . Note also that the column to column changes are more significant for the second 
root - this indicates that the "spectral fluctuations" for this root are high due to the orbit 
metamorphosis that takes place at E = V, inside of the integration interval En-i < En < En- 
This also indicates that more periodic orbit expansion terms are needed to get the more 
precise position of S2 with respect to V. 
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root 


m = 


m = 


m = 2 


m = 3 


exact 


Sl 


2.4354 


2.6198 


2.6173 


2.605 


2.5958787201295728 


S2 


6.1601 


5.5789 


5.2366 


5.1434 


4.9455316914381690 


Sl7 


53.4071 


53.405 


53.404 


53.406 


53.403119615030526 



TABLE I: The results of the periodic orbit expansion approximations to the exact values of s^. 
The first m = column gives the Weyl's Sn estimate for the roots of the spectral equation, 
and the following columns show the corrections to it due to the m = 1, m = 2 and m = 3 code 
length orbits. Note that the third periodic orbit expansion approximation to S2 comes out higher 
than the critical value, ^2 > Scrit = 5, whereas the actual root S2 is smaller than Scrit- Hence 
more periodic orbit expansion terms are needed to capture the correct qualitative behavior of 52- 

After the allowed values Sn for the action are found, the quantum levels of the momentum 
can be obtained by inverting the relationship (j36|) . 

f ^ lis <V 
k=\ 2^ ' ^' (39) 

[s if s„ < 

According to (jH) one has from ()39|): 

h = 2.605, k2 = 5.0020, = 26.4689. 

(40) 

This demonstrates how the periodic orbit expansion quantization rule (|31|) gives the explicit 
solution to our spectral problem. 



IX. A SIMPLIFICATION - SCALING POTENTIAL 

It is possible to simplify significantly the result by assuming that the hight of the potential 
step V also grows (scales) with energy, V = XE. This assumption is not as artificial as it 
may seem. Scaling (not necessarily in the V = XE form) is a common phenomenon that 



occurs in a variety of familiar physical systems [2^|23[, such as, e.g. the Hydrogen atom in 
a uniform magnetic field (also a nonintegrable system). Basically, the scaling eliminates the 
geometrical restructuring of the trajectories - i.e. in our case the transitions from classically 
allowed to "ghost" to orbits. 
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As soon as a certain value of the scaling coefficient A has been chosen, the orbits will 
either always hover above the potential step (A < 1) or sit below it (A > 1). The assumption 
V = XE also implies that the reflection coefficient Vgc (and therefore the weight factors (|15|)) 
do not depend on energy, 

r,, = 1^, (41) 

where P = Vl — A. In addition, the action lengths of the two parts of the potential (jH)) 
become simply proportional to k, Si = Lik and 5*2 = L2/3/c. The spectral equation for A < 1 
is now 

sm{QQk) = rscsin{uk) (42) 

where Qq = Li + PL2 and u = Li — [3L2 are constant frequencies, > and the roots of 
(jl^ are separated from one another by a periodic sequence of separators 

4.= j^(« + i). (43) 

For a given A, the coefficients r^c, VLq, uj and the separators do not depend on energy 
anymore. The spectral equation (jl^ doesn't change its functional form, and the tunneling 
(ghost) trajectories never appear in the system. This clearly illustrates the idea of the 
"structural freezing" of the dynamics due to the scaling. 

With so many characteristics of the system becoming constant, the evaluation of the /c„ 



greatly simplifies and the integration in can be carried out explicitly 

Let us assume that energy scales above the potential step, A < 1. The exact eigenvalues 
of the momentum in this case are 

where Sp = Sp/ k, ujp = nSp/So, and the weight factors Ap are given by (fT3j) in which r = 
This expression clarifies the overall structure of the periodic orbit expansions for 



= - ^ E io E ^ (^) ^^^vn) , (44) 



X. DISCUSSION 



We have studied an elementary example of a classically nonintegrable (quantum stochas- 
tic) system in the potential Despite the simplicity of the setup, the classical dynamics of 
the particle in the potential ® turns out to be extremely complicated, and this complexity 
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FIG. 9: The familiar schematic picture of the HeUum atom (Fig. 7a), which was also used as 
the physical model in 1920s, is analogous to using only Newtonian trajectories in the potential 
(jHl), shown in Fig. 7b. This scheme completely ignores the actual classical dynamical complexity. 
Compare the Fig. 7(b) to Fig. 3 and Fig. 5. 



is manifested in its spectral properties in quantum regime. Surprisingly, in studying this 
most elementary example, one runs into essentially all the dynamical and physical effects 
(integrability versus nonintegrability, the exponential proliferation of the periodic orbits in 
the nonintegrable case, periodic orbit expansions, the use of symbolic dynamics, phase space 
metamorphosis, tunneling, ray splitting, etc.) that appear in semiclassical analysis of more 
realistic physical systems. Usually these phenomena are investigated via a rather involved 
mathematical apparatus, whereas in our example they naturally come into play and can be 
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analyzed by elementary means. Due to its illustrative simplicity, this and similar 
systems can be thought of as the "Harmonic oscillators" of quantum chaos. On the other 
hand, this problem is rich enough to illustrate the essence of the difficulties associated with 
the semiclassical quantization of chaotic dynamical systems, such as Helium atom. 

As mentioned above, the early attempts to quantize the Helium atom failed because 
the qualitative difference in the dynamical complexity between Hydrogen and Helium were 
overlooked. The attempts to quantize a chaotic system within the framework of Bohr- 
Sommerfeld or EBK quantization theory using only a few integrable-like trajectories can 
be compared to considering just the two "naive" (Newtonian) classical trajectories in the 
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potential (L) for E < V and (LR) for E > V. From the structure of the exact result 
(jnH) it is clear that such consideration would produce only very approximate results, far 
from the real complexity of the problem. 

Formula (j31|) represents the "modification" of the EBK quantization condition (@)) men- 
tioned in the Van Vleck citation above, for the simple case of the potential (jHI), in an explicit 
and self-contained form. It would be natural to expect that obtaining the semiclassical 
spectrum in the form En = ■■■ for more complicated systems such as Helium atom should be 
a much more difficult task. However, this result creates an interesting precedent that may 
indicate new directions in the semiclassical quantization theory and related fields. 
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